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^ . A breathing soliton-like structure in dispersion-managed (DM) optical fiber sys- 

< 



tern is studied. It is proven that for negative average dispersion the breathing sohton 
is forbidden provided that a modulus of average dispersion exceed a threshold which 



depends on the soliton amplitude. 
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> : Propagation of optical pul=e in nonlinear media with varying di.per.ion is both 

, fundamentals and important applied proble because a dispersion managed (DM) sys- 
00 , 

O , tem, which is a system with periodic dispersion variation along an optical fiber, is one of the 
O _ 

O I most prospective candidate for ultrafast high-bit-rate optical communication lines. Lossless 
propagation of optical pulse in DM fiber is described by a nonlinear Schrodinger equation 

d 

j> ! (NLS) with periodically varying dispersion d{z): 

X. 

c3 ■ iu^ + d{z)utt + \u\'^u = 0, (1) 

where u is the envelope of the optical pulse, z is the propagation distance and all quantities 
are made dimensionless. Consider a two-step periodic dispersion map: d{z) = do + d{z), 
where d{z) = di for < z + nL < Li and d{z) = d2 for Li < z + nL < Li + L2, do is 
the path-averaged dispersion, di,d2 are the amplitudes of dispersion variation subject to a 
condition diLi + d2L2 = 0, L = Li + L2 is a dispersion compensation period and n is an 
arbitrary integer number. Eq. (|I|) also describes pulse propagation in a fiber with losses 
compensated for by periodically placed amplifiers if the distance between amplifiers is much 



less than L. 

In a linear regime, in which the nonlinear term in Eq. (|I]) is negligible, the periodical 
variation of dispersion is a way to overcome pulse broadening due to the chromatic dispersion 
provided that the residual dispersion do is small enough. However in real optical fiber the 
nonlinear term in is important because the optical pulse amplitude should be big enough 
to get high signal/noise ratio. One of the fascinating feature of DM system is the numerical 
observation of a space-breathing soliton-like structure, which is called DM soliton, for both 
positive and negative residual dispersion d^. This observation is in sharp contrast with the 
system described by NLS with the constant dispersion where stable soliton propagation is 
possible only for the positive dispersion0 because the nonlinearity can continuously com- 
pensate the positive dispersion only. In DM soliton the balance between the nonlinearity 
and dispersion is achieved on average over the dispersion period L what lift a restriction of 
the positive dispersion sign. Nevertheless it was never proven that DM soliton really exists 
because there is a possibility that this is rather a long-lived quasi-stable breathing pulse 
which decays on a long distance z. It is shown here that for negative d^ DM soliton can exist 
only if |(io| is small enough to allow nonlinear compensation of pulse broadening due to the 
dispersion over distance L. 

Eq. (P can be written in the Hamiltonian form lUz = where the Hamiltonian 



H 



d{z)\u,\^- 



\u 



dt, (2) 



is an integral of motion on each interval of a constant dispersion d{z) = const. Eq. (|l]) 
is reduced to usual NLS on such intervals. At points z = nL and z = nL + Li, where 
n is an arbitrary integer number, the Hamiltonian experiences jumps due to jumps of the 
dispersion although the value of m is a continuous function of z in these points. In contrast to 
the Hamiltonian the time- averaged optical power N = J \u\'^dt or number of particles in the 
quantum mechanical interpretation of NLS (in this interpretation the coordinate z means 
some "time" and actual time t has a meaning of "coordinate") is an integral of motion for 
all z. Consider z-dependence of the quantity A = J t'^\u\'^dt. A/N is the average width of a 



time-distribution of u or simply (t^) in a quantum mechanical interpretation of NLS. 
Using (|1|) and integrating by parts one gets for the first z derivative 

A, = d{z) J 2it{uu; - u*ut)dt. (3) 

In a similar way after a second differentiation by z one gets 

A,, = AdH + Ad'^X + (4) 

where X = / \ut\'^dt. It follows from Eq. (|]), which is often called virial theorem (see e.g. 
Refs.0lii), that A^ experiences finite jumps corresponding to jumps of a step-wise function 
d{z) : 

do + d2 



A.. 



=ii+o do + di 
do + di 



A, 



2=Li-0 



z=L-0 



(5) 



z=L+0 do + d2 

Set X{z) = Xo + 6X{z), X{0) = Xq then one can integrate Eq. (|) over intervals 
(0,Li), {L^,Ly. 



A. 
A. 

where 



2=Li-0 



z=L-0 



A, 



A, 



+ 4 / [(do + d,)H^ + {do + d,yX]dz 

z=0+0 Jo 

L 



2=Li+0 



4 f [{do + d2)H2 + {do + d2fX]dz, 



(6) 



Hi = {do + di)Xo - Fo, 
= {do + c^2)Xo - Fo - {di - c/2)5X 



2;=Li 



(7) 



are the Hamiltonian values on intervals (0, Li), {Li,L) respectively, Y{z) = 

I |4 

/ ^dt, Yq = Y{0). Here the conservation of Hi on interval (0, Li) is used in derivation of 
expression for H2. 

The DM soliton solution of Eq. (0) (see Ref.lll) is given by u = u{z, t) exp{zkz), where k is 
an arbitrary real constant and u{z + L, t) = u{z, t) is a periodic function of z, u{z, t)\^^^_^^ — >• 
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0. Thus for DM soliton A, 
into the form: 



=L+0 



This condition can be cast via Eq. ©J (|3) 



..J + 



f \do + dif6Xdz+ f {do + d2f6Xdz = 0. (8) 

Next step is to consider 5X{z) dependence. Using (HD and integrating by parts one can get 

X. = Aj (PtRtR'dt, (9) 

where u = Re^'^, (p and R are real, R > 0. Consider an upper bound of Xz which is given 
by a chain of inequahties: 

4 J (PtRtR^dt < 4 max(i?2) J \(l)^RtR\dt < 4X^/2 AtVs^ ^^q^ 
where the following inequalities are used: 

2(t>tRtR < {<PtR? + Rt, 

max(i?2) < 1* \{RWdt' < J \{R%\dt < 2 J R\Rt\dt < 2N^'^X^'^ (11) 

(in last expression the Cauchy-Schwarz inequality is also used). Eq. (^ and (|T0[) can be 
integrated by z and give (it is assumed below that 2Xq^'^ N^^'^max{Li, L2) < 1): 



^ ^ . 2 (12) 

(1 - 2Xo/'iVV2^)2 



In a similar way using inequality X^ > —4 / \(f)tRt\R^dt following from (^ one can get the 
lower bond of X{z) : 

X > ^° ^ (13) 

For DM soliton X{L) = Xq and thus it is more convenient to use for Li < z < L similar 
inequahties: 



Eqs. (D, (IID, (|T3D, (HD result in inequality: 



\2dQXo — lo| ^ 



oil— d2 1-^2^0 



L 



\do+di\L 



;i-2XyVl/2Li)' 
(d(l+d2)^L| 



(15) 



Eq. (|T5D is the main result of the present paper. Eq. (|I5|) is a consequence of initial 
assumption that DM soliton exist for given parameters Li,L2,dQ,di,d2 and integral values 
Xq,Yq,N which depend on u\^^q only. Thus DM soliton can exist only if this inequality is 
fulfilled. 

Note that if one assumes uniqueness of DM soliton solution for given k and soliton width 



then, as shown in Ref. I 



\u\ 



2 = 



\u\ 



. In such the term 6X 



in Eq. 



vanishes and instead of (|T^) one can get a stricter inequality. Here however that possibility 
is disregarded for the sake of generality. 

To clarify physical consequences of Eq. (|I5|) consider the optical pulse with a typi- 
cal amplitude p and a typical time-width to- Then ~ bP^o, -^o ~ bP/^o and thus 
Xq^^A^^/^L ~ L/Zni, where Z„,/ = is a characteristic nonlinear length. In a typical 

experimental condition a nonlinearity is small: L/Z^i ^ 1 and denominators in (|15D can be 
series expanded thus giving 



|2(ioXo — lo| ^ 



2\di - d2\LiL2 + \do + dilLj + 



{do 



d2?Ll 



1^0 + di 



0{ 



diL\ 



;i6) 



Provided that d^ is negative both terms in left-hand side of (|T^) have the same sign and 
thus right-hand side should be greater or equal to 2|(io|Xo + Fq- Assuming d\ ^ \do\ one can 



get from (|T^) the following estimate (Iq ~ ^o/-Z^; 



2Mo| _^ ^ < 4Lirfi ^^ ^ Li 
toZni Zli Zlito L 



(17) 



Consider a strong dispersion management limit Zdisp/L <^ 1, where Zdisp = ^diLi/tl is 
a typical dispersion length. This limit implies that an optical pulse experiences strong 
oscillation on each period L due to dispersion. Then ( [T7| ) reduces to 



do 
di 



< 




(18) 



i.e. a nonlinearity (amplitude of the optical pulse) should be strong enough to allow DM 
soliton solution existence for a given negative do. 

Eq. (|15D gives a necessary condition for DM soliton existence but not sufficient. In 
other words violation of the inequality ([T5| ) means that DM soliton is forbidden. Of course 
it would be interesting to find to what extent this necessary existence condition is close 
to sufficient one. In general this could be done only if one found DM soliton analytically. 
Here one can only mention that there is a qualitative correspondence between threshold 
of DM soliton existence following from the analytical condition (^) and from a numerical 
investigation of DM soliton. Namely the maximal value of {d^l [do < 0) for which DM 
soliton exist grows with increase of the dispersion map strength L/Zdisp according to both 
numerics (see e.g.0'0) and the analytical condition (|16|) . It also follows from ([TsD that for 
asymmetric dispersion map Li ^ L2 maximal possible value of \do\ grows as Li increase (for 
fixed L, Znl, di) in correspondence with Fig. 3 of Ref.lli. 

Eq. (|15D has also a clear physical meaning in another limit ^ ^ ^^isp ^ L and 
Znl ^ L in which ([T5| ) reduces to: 



Equality 2doXo = Yq exactly corresponds to one-soliton solution of NLS with dispersion do 
see Ref.0) where the dispersion do and the nonlinearity continuously balance each other. 
Thus in the limit Z^isp ^ L, which is called a weak dispersion limit, we recover usual 
NLS describing a path-averaged (over space period L) DM soliton dynamics provided do is 
large enough. A weak dispersion management limit was studied earlieifl'&H. Note that an 
additional condition ^ ^ allows the amplitude di of the dispersion variation still to be 



much higher than do because one assumes L <^ Zni- 

In conclusion the necessary analytical condition (|15]) of DM soliton existence is estab- 
lished. From a physical point of view this condition means that DM soliton solution can 



{2doXo - Yo)/Yo = 0{L/Zmsp) < 1- 



(19) 
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exist only if the nonlinearity is strong enough to compensate the pulse broadening due to 



the negative value of the average dispersion cIq. Note that estimates in Eqs. (|T6|) — (|T9[) 
are only given here for a physical interpretation of the analytical condition (piSl). So far 
DM soliton solution was obtained numericallyi'30, by variational and other perturbative 
approache These results are in agreement with the condition (|T3p. But analytical 

proof of DM soliton existence in the parameter region satisfying the condition (^), i.e. the 
sufficient existence condition, is still an open question. 
The author thanks I.R. Gabitov for helpful discussions. 
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